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o , 

^ A new method describing nuclear rotational motion microscopically is proposed. We 

CNJ , extract the rotational Hamiltonian by introducing the intrinsic pair modes which commute 

with the rotational mode. Thereby the rotational mode is not treated as zero energy mode 
in contrast with the conventional RPA formalism so that we circumvent the difficulty related 
with infrared divergence. The wave function is constructed by angular momentum projection 
on each intrinsic state. Without numerical integration for projection we calculate the matrix 
' elements analytically under a certain approximation. The numerical calculations are carried 

out to illustrate the applicability of our method and they show that our method works well. 
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§1. Introduction 



The study of rotating nuclei is one of the main topics of the nuclear physics. In 
order to describe nuclear rotational motion microscopically, the mean field approxi- 
mation is often used, such as the cranked-Hartree-Fock-Bogoliubov (CHFB) theory^\ 
\ Although this approach succeeded in explaining the various phenomena of rotating 

nuclei, it has a serious defect; the method inevitably accompanies the spontaneous 
symmetry breaking and thus the angular momentum is conserved only on the average. 
CJ \ The resultant wave functions are no more the eigenstates of angular momentum. 

Many works toward the restoration of broken symmetry have been made. One of 
the useful approaches is the projection method'^^~^\ that is the angular momentum 
projection (AMP) in this case. AMP is a mathematically rigorous method. How- 
^ ■ ever it needs much of a numerical effort. Furthermore the conventional treatment 

of AMP loses the simple picture of intrinsic excitation. Therefore the projection 
method has been applied in approximate version . Owing to the recent devel- 
opment of computer technology, the calculations with full projection, that is the 
projected shell model (PSM) calculations (for a review, see Ref. 9), has been car- 
ried out although they are still limited. The application to axially deformed nuclei 
are carried out ^'^^ " . Furthermore the calculation including triaxiality have been 
made " ^'^^ in transitional region. In these works the Nilsson-plus-BCS state is 
adopted as the intrinsic wave function, that means the self-consistent treatment of 
the coupling between the rotational mode and intrinsic mode is not fully considered. 
The application to the CHFB state, in which the axial symmetry is no more kept, is 
also made^^^ -20)^ ^j-^g effect of the intrinsic excitation is not took into account. 
The random-phase approximation (RPA) is another useful approach. This method 
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separates the rotational mode as zero-energy mode from the intrinsic mode and re- 
stores the broken symmetry within its approximation. However this method also has 
a shortcoming ^-"^^ ; it is valid for only the small amplitude oscillation, while the zero 
energy mode ranges over the large amplitude because of no restoring force. This 
discrepancy, for example, appears as the divergence of the norm of RPA ground 
state. Although the attempts to overcome this difficulty have been made, the de- 
scription of finite rotation starting from the small amplitude approximation has its 
own difficulty " . 

In this work we propose a new method to get rid of the difficulty of RPA and 
to give an approximate treatment for AMP at the same time. The essential point is 
that we do not treat the angular momentum as an eigenmode although the intrinsic 
mode is decoupled with it. The wave function having the definite angular momentum 
is formally given by the angular momentum projection. However we give an analytic 
expression for the matrix elements between these states without integration. For 
simplicity, in the present work, the nucleus is assumed to be axially deformed. 

In §2, the outline of our method is shown together with a brief review of RPA. 
§3 is devoted to the formulation of our method in which the rotational Hamiltonian 
is extracted by introducing the intrinsic operator. The wave functions are obtained 
in the subspace spanned by the multi-intrinsic pair base. In §4, we derive the ana- 
lytic expression for the moment of inertia and B{E2) transition probability under a 
certain approximation. The numerical calculation is made to check the validity of 
our method and the results are compared with those of RPA and AMP in §5. Finally 
the summary is given in §6. 



§2. Basic Idea 



We begin this section with a brief review of RPA. In axially symmetric system, 
Jz has the good quantum number. The RPA vibrational phonon operators Ojj, On 
are derived from the equation 



RPA 



EnOl 



(2-1) 



The subscription RPA denotes that the only leading terms are considered. In the 
solution of this equation, however, the zero energy mode J± appears 



H,J± 



RPA 



0, 



(2-2) 



where J± = Jx ^ iJy Eq.(2-2) means that the symmetry broken by the mean field 
approximation is restored. Hrpa is completely divided into the rotational part and 
vibrational part within its approximation. 

Hrpa = Erpa + + ^ EnOlOn, (2-3) 

where Erpa = {Orpa \H \Orpa ) and \Orpa ) is RPA vacuum defined by 

On\0RPA) = 0, (2-4) 

J±\Orpa)=0. (2-5) 
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However this definition leads to the well-known diflaculty; the norm {Orpa\Orpa) 
diverges infinitely. It stems from the fact the rotational mode is treated as an 
eigenmode. 

Our approach is proposed in order to avoid this difficulty. We also separate 
the rotational mode and intrinsic mode but we never treat the rotational mode as 
an eigenmode. For the first step, the intrinsic pair operator is defined as the 
commutable mode with J± in lowest order, 

J±,Xl\=0. (2-6) 

The second step is to represent the Hamiltonian in terms of J±,Xj^,Xn and to 
decompose it to three parts, 

H = Hrot + Hintr + Hcoupl + (const.), (2-7) 

where 

Hrot = Hrot{J-J+), (2-8) 

Hintr = HintriX]^, Xn), (2-9) 

Hcoupl = Hcoupl{J±,Xl^,Xn)- (2-10) 

For the last step, H is diagonalized in the subspace spanned by the projected multi- 
intrinsic pair base 

{\IM)} = {PU |0) ^PLkXI \Q),PUxixl |0) , ...}, (2-11) 

where Pmk is the angular momentum projection operator. 

In the usual projection method like PSM, H is diagonalized in the projected 
Tamm-Dancoff base 

{\IM)} = {PU |0) ,P[,K (a^a^. '^^^ ("^"0. ' (^-^^^ 

It should be noted that {a^a^).,^, (a^a^)i/ include both the rotational and intrinsic 
modes, while the rotational mode is already eliminated in X^ in lowest order. This 
leads to an advantage of our approach that the analytic form of matrix is approxi- 
mately obtained without integral calculation. (See §3.) 

§3. Formulation 

3.1. Intrinsic Pair Operator 

In axially symmetric system, any operator is classified by K quantum number. 
In this case, J± corresponds to the rotational mode. It is written in the quasi-particle 
representation as 

J±=Y1 {-^fapaial + Jfo^paf^aa} + J±ai3aia/s, (3-1) 

a</3 a/3 
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where we follow the notation of Ref. 1). The collective pair operator X^. is defined 



as 



= E ^IS^Ul (3-2) 



1/3 "-a"-/?! 
a</3 

J± = ao {Xi + X^) + E 4 (3-3) 

a/3 

~ " 1/2 

where oq = (0 | J_ J+ |0) . Since oq is large enough in deformed nuclei 

4-ao(xi+X^) (3-4) 

The intrinsic pair operator is defined so as it may satisfy (2-6) up to the leading 
order, that means 

a<l3 

E = 0- (3-6) 

a</3 

Prom now on, we use the following notation rule; the subscription 'it' or 'c' 
denotes the collective mode and 'n' denotes the intrinsic one and '/u, i/, denotes 
both. If necessary, K quantum number subscription is added, as e.g. X^j^. The 

amplitudes tp'^p satisfy the orthogonality and completeness 

E ^S^S = V, (3-7) 

a</3 

E = <^-a' V (a < /?). (3-8) 

They can be simultaneously determined from the eigenvalue equation for the Stability 
matrix 



where 



E E ^aT('^"ft"'/3'+'^a/3;^)V'S' =^Ai'' (3-9) 
a</3 a'</3' 



^a/3;a'/3' = (0 |[a/3aa, 4,aJ,]] |0), (3-10) 

'Sa/3;a'/3' = -(0|[a/}O«, [i/, Oa'a/j/]] |0), (3-11) 



and a' (3' denotes the time-reversal index of a'P' . 

3.2. ExtracfAon of the rotational Hamiltonian 

The Hamiltonian with two-body interaction is written in terms of pair operators 
up to leading order (Appendix A), 

H ~ {H) + Y,^nXiXn + ^ {Xlxl + X^X^ + 2Xtx,)] , (3-12) 

n (j,v 
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where Vy(^ni,^ is a newly defined matrix 

= E E ^ifHt%a','^%. (3-13) 

af3 a'P' 

By classifying (3-12) with respect to each mode, H can be decomposed into the 
rotational part Hrot which can be written only by the collective pairs, the intrinsic 
part Hintr only by the intrinsic pairs, and their coupling part Hcoupl 

Hrot = ^%^J_ J+, (3-14) 

Hintr = E ^nXlXn + E Vv(nn') (^^4' + ^n'^n + 2XtX„,) , (3-15) 

n n,n' 

Hcoupl = 2 E ^v{nc) {x^xl + XcXji + Xl^Xc + xlXfi^ . (3-16) 

n,c 

In Eq.(3-14) we used (3-4). 

3.3. Wave function 

We span the intrinsic space 

\ki) = Xl...xi^\(}) , (3-17) 
and make the wave function in the following base 

{\IM)} = {Pi,K\k)} (3-18) 

where 

^^K = -^l daJ^dPsinPj^ djD'MKic^,/3,^)e-"''^e-'^'ye-'^'^. 

(3-19) 

Because J± change into c-numbers when they are operated to these states, all we 
have to do is calculate the norm matrix of this base. 

Here we note that K > intrinsic pairs have to be prepared to compose the 
orthogonal base because 

(0|X,^P^^Xt^|0) = {0\X,,^KPLKKXiK\0) 

= {OIX^.kPLk-kXI.k |0>. (3-20) 

This property stems from the time-reversal symmetry of the HFB wave function 
in the axially deformed system. Moreover, through this study, we neglect the non- 
diagonal elements of the norm matrix 

{0\PLXIJ0), {0\Pi,KXixl\0), (0|X„Pi,^4,4„|0),... (3-21) 

compared with (0 |ioo |0). This is because Pmx includes only the rotational mode 
which is orthogonal to the intrinsic mode in lowest order. Under this approximation, 
we can calculate the matrix without integration. 
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§4. Physical Quantities 

4.1. Moment of Inertia 

4.1.1. Peierls-Yoccoz' Moment of Inertia 

The Peierls-Yoccoz' moment of inertia^^^ is defined 



Ei = {IM\H\IM) 

2JpY 



E.+ '4P^, («) 



where \IM) is determined as following 

|7M)=<Pi,o|0) (4-2) 

where the normalization constant A/'q = (OlPoolf*) ^^'^ ■ Comparing (4-1) with 
(3-14) — (3-16) under the approximation mentioned above, the Peierls-Yoccoz' mo- 
ment of inertia can be derived analytically 

Jpy = TTF^- (4-3) 

4Fi;(+_) 

4.1.2. Thoulcss-Valatin Moment of Inertia 

As a rule the Peierls-Yoccoz' moment of inertia tends to underestimate the ex- 
perimental value and thus we extend the physical space by adding K = 1 one 
phonon states perturbatively 

\IM)^NiPU\^)+r^ K+^PLxIA^)^ (4-4) 

nK=+i 

where 



K^ = {0\Xn,PkA\Q) ' 

The perturbative energy shift in second order is given by 



AE = - (4.5) 



Here we must mention that the definitions of uJn and Xl^^ are slightly different from 
those of (3-9) because the formers are defined as the Tamm-Dancoff eigenvalues and 
eigenvectors in the intrinsic space, i.e. they are derived by diagonalization of A after 
removing the zero energy modes. 

In addition to (3-21), we also neglect the terms 

(0|X„^,P/i<;JO). {n^n') (4-6) 

compared with (0 |Poo |0)- can be derived analytically 

AEc^-M{I + l)af Y <X\n+i-). (4-7) 

'nK=+i 
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Here we used the approximate relation (see Appendix B), 

(0|X„^,P/i^l+jO) K(+n_i) _ 1 



(Ol^o'olO) 



2K(ri+i-) 2- 



(4-8) 



which can be naturally derived from our approximation. Thereby the moment of 
inertia is newly defined 



(4-9) 



where J^eff corresponds to the Thouless-Valatin moment of inertia JtVi which co- 
incides with Jq in (2-3). It is calculated from 



^0 = 2 E V'it(^-^)aL'/J'V'S'- 



(4-10) 



4.2. B(E2) Transition Probability 

The B{E2) transition probability in the ground band also derived in the same 
framework. In the space \IM) = J^qPmo |0), it can be written 

B{E2;I + 2^ I) 



KrI+2 2 

-ijT E (/-A;2A|/ + 20)5a(/) 

•'^O A=-2 



, (4-11) 



where g\{I) is given by 



9x{I) = (-1 



,a (0|^oaQ2a|0) 
(OjPo'olO) 



(4-12) 



g\{I) and Mq^'^ /Mq can be given the analytic form in our approximation. (See 
Appendix C.) 



9x{I) - < 



Q20) (for A = 0) 

-e±ao V^-f + l) (for A = ±l) 
(for A = ±2) 



(4-13) 



Here is the overlap coefficient 

W = (0|^MKQ2;f|0) (4-14) 

and 

(0|Pjo|0) / 2/ + I Ea(-^+2A;2 -A|J0) 9^+2) 
[Xi J (0|Po'o+'|0) V2/ + 5 Ex{I\;2 -X\I+20) g^^i) ' 

(4-15) 

Prom (4-15) and the normalization condition J2i (0 1-^^ |0) = 1) the angular momen- 
tum distributions of the HFB wave function |0) can be analytically evaluated. 
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§5. Numerical Calculation 

We carry out the numerical calculation applying our method to Er isotopes and 
compare the results with AMP and RPA. As the purpose is to check the validity of our 
method, we adopt a rather simple effective interaction; the Pairing-plus- Quadrupole 
interaction. The detail correspondence with experimental data by using the more 
realistic interaction is left in the future study. 

The adopted Hamiltonian is 

1 ^ 

H = Ho--xYl Q2XQ2X - GP^P, (5-1) 

^ A=-2 

where the chemical potentials are already absorbed in the spherical single particle 
energy part Hq. The single particle basis employed is A'^ = 4,5 and 6 for neutrons (64 
levels) and N = 3,4 and 5 for protons (46 levels). The exchange terms are neglected. 
The force parameters are below; for the Q-Q interaction 

Xnn = 3.768 X 10-2 (MeV), (5-2) 
Xpp = 4.796 X 10-2 (MeV), (5-3) 

Xnp — Xpn — V XnnXpp > (^"^) 



where Q is scaled to be dimensionless by the oscillator length 60 = \ tttt- For the 
pairing interaction we use 

N - Z 

G = (22.36 ^ 14.59 )M (MeV), (5-5) 

where - for neutron and + for proton. 

In table I and II, the calculated results for the moment of inertia JpyiJt.v. are 
shown. They indicate that our method gives a good approximation for AMP and 
RPA. 



Table I. Comparison of Jp.y.^ • Table II. Comparison of Jt.v.^ ■ 

MeV"^) obtained from AMP with MeV"^) obtained from RPA with 

the one from our method the one from our method 



le^Er ^«^Er ^*'*'Er 
AMP 20.6 22.7 24.5 26.8 
ours 22.4 24.5 26.1 28.6 



i«^Er ^''^Er ^««Er 
RPA 29.8 33.1 35.0 38.0 
ours 31.3 34.9 36.8 39.8 



The small discrepancies among them are supposed to be the consequence of 
our approximation (3-21). In order to check the accuracy of the diagonal element 
(0 l^n+i A^i^i+i |0 ) /(O |Pjo |0 ) (4-8), we calculated it by the exact AMP. The results 
are shown in Table III. Here we showed the data for the lowest energy intrinsic mode. 
(We checked the results for the other mode are nearly same.) 
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The B(E2) transition probability and the angular momentum distributions (0 |0 ) 
for ^^^Er are shown in Table IV and Fig.l. In this calculation, we used the isoscalar 
effective charge Ceff = 0.3e. Our method can well reproduce the result of AMP. 
To look into the origin of good agreements, we calculated gx{I) by the exact AMP. 
The results are given in Table V. We find that the main part of B{E2) comes from 
gx=oiI), which does not depend on the detail of the projection because it is mainly 
determined by < Q20 >■ The same reason holds for the good agreement of (0 |Pqo |0) 
with the exact one. 



Table III. Comparison of Eq.(4-8) 
obtained from AMP with the 
one from our method 



B(E2;I+2 ^ I) (e^ • 6^) of 
Er obtained form exact AMP 
and our method. 



Table IV 

166 



I 


AMP 


approx. 


I 


AMP 


ours 


2 


0.49184 







1.55 


1.59 


4 


0.49313 




2 


2.22 


2.26 


6 


0.49515 


0.5 


4 


2.44 


2.48 


8 


0.49788 




6 


2.56 


2.59 


10 


0.50129 




8 


2.63 


2.64 


12 


0.50536 




10 


2.68 


2.67 








12 


2.71 


2.68 




12 16 20 24 28 32 36 40 
Angular Momentum I 



Fig. 1. Angular momentum distribution (0 |0) for ^®®Er 



Table V. g\{I) values obtained from exact AMP and our approximation 





9\= 


=o(/) 


9\= 




9x 


MI) 


I 


exact 


approx. 


exact 


approx. 


exact 


approx. 





77.6 















2 


77.6 




1.62 


1.61 


0.014 




4 


77.7 




2.96 


2.94 


0.053 




6 


77.9 


78.5 


4.29 


4.26 


0.113 





8 


78.1 




5.61 


5.58 


0.195 




10 


78.3 




6.93 


6.90 


0.297 




12 


78.6 




8.24 


8.21 


0.419 
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§6. Summary 

In this work, we proposed a new niieroscopie method for the description of nu- 
clear rotational motion that avoids the difficulty of RPA. We separate the rotational 
Hamiltonian by defining the intrinsic mode that commutes with the rotational mode 
up to leading order. The wave function with definite angular momentum is formally 
determined through angular momentum projection, though practical integral calcu- 
lation is replaced by the analytic estimation under a certain approximation. We 
calculated moment of inertia and B[E2) transition probability in the ground band 
by using this method, restricting ourselves to axially deformed system. The obtained 
results shows good agreements with the ones from AMP and RPA. 

This method has various advantages as following; (i) It keeps the clear physical 
picture of RPA treating the rotational and the intrinsic modes separately, (ii) all 
the necessary calculations are reduced to merely analytic ones and thus it is easier 
to include much of an intrinsic excitation. 

Finally we mention that this work is still preliminary. To compare with exper- 
imental data, more realistic interaction should be applied, such as the pairing-plus- 
Quadrupole-plus-Quadrupole-pairing interaction, the Gogny interaction. Moreover, 
more elaborate formulation that deals with K quantum number mixing has to be 
made for the study of high spin state. These works are now in progress. 
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Appendix A 

Phonon Representation of Hamiltonian 

The Hamiltonian with two-body interaction can be written in the quasi-particle 

space 

H={H) + H^^ + ^(^20 + h.c.) + ^H^^ + {H^^ + h.c.) + {H"^ + h.c), 

(A-1) 

where 

{H) = {0\H\0), (A-2) 

a 

a/3 a'l3' 
a/3 a'/3' 
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a/3 a'P' 

and H"^^ = in the HFB ground state. Ea is the quasi-particle energy and Bap is 
the scattering operator defined by 

Ba/3 = ala/S, (A- 7) 

Between Ba/s and the pair operators, there exist the following commutation relations 

[X^,Xl]=5,^ + J2r^j'^B^f,, (A-8) 

a/3 

[B^/s,xl] = -Y,rifxt, (A-9) 

1/ 

where 

7 

is the first order of magnitude with respect to the small parameter e 

^ Vo the 

number of states). B^p can be rewritten 

4/3=-E^ir^4^-+o(e') (A-ii) 

since this form satisfies (^4-9) except for O(e^) deviation. Using the newly defined 
matrices 

= E E ^ifHl%^'P'^%, (A.12) 

Vvi,.) = E E V'?/^a°^;a'/3'^S', (A-13) 

Q/3 a'/3' 

^2/(M;a'/3') = E V'?/^a]3;a'/3', (A-14) 
a/3 



if is rewritten as 
H = {H) 

+ E I (- E + i^-(M-)) 4^- + ■^kmp) (44 + ^p^^) | + o{e% 

(A-15) 

As the Stability Matrix is (without exchange terms) 

-^a/Sia'/S' = {Ea - E/3){Saa'Spi3' - ^ap'^Pa') + -f^a|;a'/3'> (A-16) 
BaP;a'l3' = ^H^fj-a'fS'^ (A" 17) 
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(3-9) can be rewritten 

Prom (A- 15) (A- 18), the phonon representation of H (3-12) is derived. 

Appendix B 

Evaluation of the AE 



The perturbative energy shift in the second order is 

ae=-y: (B-i) 



Prom 



{0\PiMHPLiXi^, |0) = {0\HPi,Xi^, |0) (B-2) 

= (0|Po'i^iJXt^JO). (B-3) 

Calculating {B-2), {B-3) separately under our approximation, 

{B-2)c^{0\H,oupiP^iXl^,\0) 

~ 2K.(n+i -)(0|X„^,X_Po'i^i+, |0) +2y„(n_i +)(0 |X„_,X+Po^iXt^^ |0) 
= V^(/ + l) 4ao-V,(n+i -)(0|X„^,P/i^i+i |0), (B-4) 

(5-3)~(0|Po'ii^co«piAl+-JO) 

:^ ^/(/ + 1) 2ao-V,(+ n_i)(0 |Po'o |0). (B-5) 

Prom {B-A) = {B-5), 

(0|X„^,P/i^l+jO) K(+n_i) _ 1 



(0|Po^o|0) 2K(n+i-) 2' 



(B-6) 



Here we used the relation Vy{n-i +) = K(+ n_i) = T^(n+i — ) = Vy{— n+i) . As 
the result, AE can be given analytically 

AE . -/(/ + 1) E ^^^^^ (B-7) 

-8I{I + l)ao 2 E -)• (B-8) 
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Appendix C 

Evaluation of the B{E2) and (0 |0; 



B{E2) transition probability is given by 



B{E2;I+2 ^ I) = 



1 



21 + 5 



1(^+2 IIQ2II/; 



(C-1) 



where {I+2\\Q2 \\I) is the reduced matrix element derived from the Wigner-Eckart 
theorem, 

{I+2M+X\Q2x\IM) = {2I + iy^^^ {I+2M+X;2 - A |/ M ) (7+2 IIQ2 ||/>. 

(C-2) 

We have only to calculate {1+2 M \Q20\IM). In the space \IM) = MqPmq |0 ) 

{I+2M\Q2o\IM) 
lu-i+2 2 
= '-^(-l)^(/M;20|/+2 -M) J2 (/-A; 2 A |/+2 0) 5a(/), (C-3) 



where 



9xil) ^ (-1) 



A=-2 



x {0\PLQ2x\0) 
(0|Pjo|0) 



(C-4) 



Here we used the relation 

PkMQ2X= E (l'M';2\\IM){l'K';2X'\IK)Q2yP^,M'- (C-5) 



X'l'M'K' 

Finally we can rewrite the B{E2) 



B{E2;I+2 I) 



(I -A;2A|/+20) g^{I) 



(C-6) 



gx{I) and Mq'^'^/N'q can be derived approximately. Pair operator representation of 
Q2a|0) is 



[ f(Q2o)+Enen^i)|0) (for A = 0) 

Q2X\0) = 1 (e±^i+En^n^^)|0) (for A = ±l) 

, EnCnXllO) (for A = ±2) 
Prom our approximation, gx{I) can be analytically written. 

' (Q20) (for A = 0) 

9x{I)-{ -C±ao V^(-^+ 1) (for A = ±1) 

(for A = ±2) 



(C-7) 



(C-8) 
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On the other hand, AAq+V-^o^ is derived from (0 IP^^ Q20 Pmo |0) 

(0|Po'^2Q2oPi^o|0) 
= (7M;20|7+2M)(0|Po^o|0)E(-^^;2 -X\I+2 0) gx{I) (C-9) 

A 

= (I+2M;20|/M) (0 |PoV^ |0 ) ^ (1+2 A ; 2 - A |/0 ) 5a(^ + 2). (C-10) 

A 

Prom (C-9) = (C-10), (4-15) can be analytically given. 

References 

1) P. Ring and P. Schuck, The Nuclear Many-Body Problem (Springer- Verlung, Heidelberg, 
New York, 1980). 

2) K. Kara and S. Iwasaki, Nucl. Phys. A332 (1979), 61. 

3) K. Hara, S. Iwasaki, and K. Tanabe, Nucl. Phys. A332 (1979), 69. 

4) K. Hara, and S. Iwasaki, Nucl. Phys. A348 (1980), 200. 

5) K. Hara, A. Hayashi, and P. Ring, Nucl. Phys. A385 (1982), 14. 

6) K. Hara, and S. Iwasaki, Nucl. Phys. A430 (1984), 175. 

7) A. Kamlah, Z. Phys. 216 (1968), 52. 

8) Y. Nogami, and I. J. Zuckcr, Nucl. Phys. 60 (1964), 203. 

9) K. Hara, Y. Sun, Int. J. Mod. Phys. E4 (1995), 637. 

10) Y. Sun, J. L. Egido, Nucl. Phys. A580 (1994), 1. 

11) Y. Sun, M. Guidry, Phys. Rev. C52 (1995), R2844 

12) Y. Sun, K. Hara, Phys. Rev. C57 (1998), 3079. 

13) V. Velazquez, J. Hirsch, Y. Sun, M. Guidry Nucl. Phys. A653 (1999), 355 

14) G. Long, Y. Sun, Phys. Rev. C63 (2001), R021305 

15) J. Sheikh , K. Hara, Phys. Rev. Lett. 82 (1999), 3968 

16) Y. Sun, K. Hara, J. Sheikh, J.Hirsch, V. Velaquez, M. Gidry, Phys. Rev. C61 (2000), 
064323. 

17) J. Sheikh , Y. Sun, P. Palit, Phys. Lett. B507 (2001), 115. 

18) K. Enami, K. Tanabe, N. Yoshinaga, Phys. Rev. C59 (1998), 135. 

19) K. Enami, K. Tanabe, N. Yoshinaga, Phys. Rev. C59 (1999), 2494. 

20) K. Enami, K. Tanabe, N. Yoshinaga, Phys. Rev. C61 (1999), 027301. 

21) E. R. Marshalck, and ,J. Wcneser, Ann. Phys. 53 (1969), 569. 

22) E. R. Marshalek, Nucl. Phys. A266 (1976), 317. 

23) E. R. Marshalek, Nucl. Phys. A275 (1976), 416. 

24) M. J. A. de Vogt, J. Dudek Z. Szymaiiski, Rev. Mod. Phys. 55 (1983), 949 

25) J. Yoccoz, Proc. Phys. Soc. (London) A70 (1957), 388 



